Search for: All records

ABSTRACT In recent years, cosmological hydrodynamical simulations have proven their utility as key interpretative tools in the study of galaxy formation and evolution. In this work, we present a comparative analysis of the baryon cycle in three publicly available, leading cosmological simulation suites: EAGLE, IllustrisTNG, and SIMBA. While these simulations broadly agree in terms of their predictions for the stellar mass content and star formation rates of galaxies at $$z\approx 0$$, they achieve this result for markedly different reasons. In EAGLE and SIMBA, we demonstrate that at low halo masses ($$M_{\rm 200c}\lesssim 10^{11.5}\, \mathrm{M}_{\odot }$$), stellar feedback (SF)-driven outflows can reach far beyond the scale of the halo, extending up to $$2\!-\!3\times R_{\rm 200c}$$. In contrast, in TNG, SF-driven outflows, while stronger at the scale of the interstellar medium, recycle within the circumgalactic medium (within $$R_{\rm 200c}$$). We find that active galactic nucleus (AGN)-driven outflows in SIMBA are notably potent, reaching several times $$R_{\rm 200c}$$ even at halo masses up to $$M_{\rm 200c}\approx 10^{13.5}\, \mathrm{M}_{\odot }$$. In both TNG and EAGLE, AGN feedback can eject gas beyond $$R_{\rm 200c}$$ at this mass scale, but seldom beyond $$2\!-\!3\times R_{\rm 200c}$$. We find that the scale of feedback-driven outflows can be directly linked with the prevention of cosmological inflow, as well as the total baryon fraction of haloes within $$R_{\rm 200c}$$. This work lays the foundation to develop targeted observational tests that can discriminate between feedback scenarios, and inform subgrid feedback models in the next generation of simulations. 

ABSTRACT The scatter about the mass-metallicity relation (MZR) has a correlation with the star formation rate (SFR) of galaxies. The lack of evidence of evolution in correlated scatter at z ≲ 2.5 leads many to refer to the relationship between mass, metallicity, and SFR as the Fundamental Metallicity Relation (FMR). Yet, recent high-redshift (z > 3) JWST observations have challenged the fundamental (i.e. redshift-invariant) nature of the FMR. In this work, we show that the cosmological simulations Illustris, IllustrisTNG, and Evolution and Assembly of GaLaxies and their Environment (EAGLE) all predict MZRs that exhibit scatter with a secondary dependence on SFR up to z = 8. We introduce the concept of a ‘strong’ FMR, where the strength of correlated scatter does not evolve with time, and a ‘weak’ FMR, where there is some time evolution. We find that each simulation analysed has a statistically significant weak FMR – there is non-negligible evolution in the strength of the correlation with SFR. Furthermore, we show that the scatter is reduced an additional ∼10–40 per cent at z ≳ 3 when using a weak FMR, compared to assuming a strong FMR. These results highlight the importance of avoiding coarse redshift binning when assessing the FMR.